Stirring by swimming bodies 
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We consider the stirring of an inviscid fluid caused by the locomotion of bodies through it. The 
swimmers are approximated by non-interacting cylinders or spheres moving steadily along straight 
lines. We find the displacement of fluid particles caused by the nearby passage of a swimmer as 
a function of an impact parameter. We use this to compute the effective diffusion coefficient from 
the random walk of a fluid particle under the influence of a distribution of swimming bodies. We 
compare with the results of simulations. For typical sizes, densities and swimming velocities of 
schools of krill, the effective diffusivity in this model is five times the thermal diffusivity. However, 
we estimate that viscosity increases this value by two orders of magnitude. 



Munk [T] was the first to ask whether biology has an 
important impact on mixing in the oceans. Since mix- 
ing affects the global circulation and stratification of the 
oceans, it is of great interest to physical oceanographers 
to settle this question. Dewar et al. [5] have proposed 
that the mechanical energy delivered by the swimming 
motions of the marine biosphere could amount to almost 
10 12 W, a figure comparable to the energy delivered by 
the winds and tides. This suggests a biological origin of 
about 33% of the mixing in the oceans, an enormous fig- 
ure. By assuming that this energy is delivered to the top 
three kilometers of the oceans, they estimate an effective 
diffusivity produced by swimmers to be approximately 
.2 cm 2 / sec, or about 100 times the molecular value for 
heat. Kunze et al. [3] have measured elevated levels of 
ocean turbulence due to swimming krill, though some 
have questioned whether this turbulence can efficiently 
overturn a stratified medium [H[S]. Katija and Dabiri [BJ 
suggest that the displacement of fluid particles by swim- 
ming bodies, which viscous effects can lengthen, is more 
relevant to stirring than the scale of turbulence they pro- 
duce. This is the viewpoint adopted in this Letter. 

Huntley and Zhou [Tj considered the energy produced 
by 11 representative species of schooling animals, from 
krill to whales, and found that regardless of size the en- 
ergy input to the ocean per unit mass from the swim- 
ming of these animals was roughly a constant of order 
10~ 5 W/ kg. If the average the biomass of the oceans 
has a volume fraction c over a volume V in cubic meters, 
we can thus (assuming biological materials have the den- 
sity of water) arrive at a total energy input of 10~ 2 cV W. 
The area of the oceans is about 3.6 x 10 14 m 2 , and if we 
assume the biomass is evenly distributed to a depth D 
meters, we get an input of 3.6Dc x 10 12 W. This suggests 
that for a typical depth of one kilometer the volume den- 
sity of biomass should be of order 10 -3 -10 -4 , and thus 
that the organisms form a dilute suspension, even if their 
distribution is patchy and organized into schools. 

The object of the present paper is to determine the ef- 
fective diffusivity and the statistics of the concentration 



field of a passive scalar that can result from the fluid mo- 
tions caused by such dilute arrangements of swimming 
animals. The scalar could be for example heat, salt, or 
nutrients, and is assumed to have negligible feedback on 
the flow (at least at small scales). The animals considered 
by Dewar et al. [5] and Huntley and Zhou [7] have large 
Reynolds number, typically 10 2 -10 7 . Our focus is there- 
fore distinct from the mixing that can occur from dense, 
tightly interacting suspensions of Stokesian swimmers, as 
in [8lU3|. though it shares many common features such 
as linearity in the density of swimmers. Real swimmers 
present a wide variety of motions and wakes, but these 
have generic forms in the far field [14] . For example, the 
far potential field of a neutrally buoyant fish in steady un- 
accelerated swimming decays at least like a quadrapole. 
Following [BJ, we will model the swimmers by identical 
cylinders or spheres moving in potential flow, but these 
are merely stand-in examples. We emphasize that though 
the motivating application comes from oceanography, the 
simple model we introduce can be applied to a range of 
systems, such as the mixing caused by vortices or sus- 
pensions of solid particles. 

Dilute suspensions of swimmers. As an extremely sim- 
ple model of stirring by swimmers, we consider a fluid 
particle, called the target particle, which is influenced by 
the occasional passing of swimming bodies. We assume 
that every swimmer moves at a fixed speed U, over dis- 
tances O(L) large compared to the "range of influence." 
This range will typically be a few body lengths normal 
to the swimming path, and represents the distance from 
a target particle where interaction with the swimmer be- 
comes significant. We assume that encounters of our tar- 
get particle with swimmers are occasional, in that at each 
encounter the particle moves by a distance A (a), where 
the impact parameter a is the perpendicular distance of 
the initially unperturbed particle from the extended line 
of motion of the approaching swimmer. Since we are as- 
suming that L is large compared to the distance of signif- 
icant interaction, we may assume that by an "encounter" 
we mean that the A(a) may be computed from motion 



FIG. 1: A swimmer coming a distance a from a particle at O. 

of the swimmer along a doubly-infinite line. Each en- 
counter, the fcth say, moves the particle a distance A(a^), 
in the direction of the unit vector r^. Thus after M en- 
counters the position of a target particle initially at Xq 
is Xm given by 



XM = Xq 



M 

E 

k=l 



A(a k )r k . 



(1) 



Assuming infrequent encounters, we regard the dfe,rfe 
as independent and identically distributed random vari- 
ables. We then follow Einstein's derivation [TS] for com- 
puting the displacement, 



(M 2 ) = M(A 2 (a)) + Q2 A(a,)A(a fe ) fj ■ r k 



(2) 



where we assume spatial homogeneity to set Xq — 0, and 
the angle brackets denote ensemble averaging. The r k 
are isotropically distributed, so the second term above 
vanishes. 

To evaluate M(A 2 (a)), consider a swimmer at a large 
distance L from the target particle at the origin O (see 
Fig. [IJ. Imagine a "target disk" (or sphere in 3D) of 
radius R, with I -C R -C L, where £ is a typical length 
scale for the swimmer. The fraction of swimmers that 
will hit the target disk (sphere) is 2R/2ttL (it R 2 / 4tt L 2 
in 3D). We can use the number density and the volume 
of the 'shell' at a distance L to find that the number of 
swimmers hitting the target from a distance L is 2Rn dL 
{irR 2 n dL in 3D), where n is the number density of swim- 
mers. We now integrate from to Ut, since swimmers 
further than Ut cannot hit the target, to find the number 
of swimmers M that will hit the interaction disk (sphere) 
in time t: we find M = 2RUnt (M = TrR 2 Unt in 3D). 
The expression for the squared displacement ^ is now 



M 



p(a)A 2 (a)da 



since the largest value of the impact parameter is R. 
Here p{a) is the probability distribution of impact pa- 
rameters; since the swimmers are assumed to arrive from 
far away, this is p(a) = 1/R (p(a) = 2-Kaj-KR 2 in 3D). 
Combining these results and taking R —¥ oo, we find 



x = 



_ hUntf™ A 2 (a)da (2D); 
\2nUntJ™aA 2 (a)da (3D), 



(3) 
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FIG. 2: Displacement A(a) due to a cylinder, as a function of 
the impact parameter a. The dashed lines are the asymptotic 
forms for small and large a of Eq. (J5J) . A few typical trajecto- 
ries are overlaid: for large a the trajectory is almost circular, 
but is still ribbon-shaped. 

assuming that the integrals converge. Since the effective 
diffusivity k is defined by (|a;| 2 ) = 2dnt, with d the spatial 
dimension, Eq. ([3| can be used to determine n. 

Displacement due to a moving body. To find the dis- 
placement A (a) of a target particle due to a swimming 
body coming from infinitely far away and swimming in 
the z direction, we need to integrate the equations for 
the position of a target particle at x(t), 



x 



u(x,t) , x(— oo) = (a, 0, — oo) 



(4) 



for t ranging from — oo to oo. The impact parameter, 
a, appears as the initial vertical position of the target 
particle. This classical 'drift' problem has been treated 
in great detail by many authors; see for example |16H19j . 

As a working example of the displacement due to a 
moving body, we shall first treat the cylinder. We con- 
sider the displacement in two-dimensional flow due to the 
passage of a cylinder of radius £ moving at speed U in 
potential flow. Equation Q can be integrated numeri- 
cally for a given a to obtain A (a), plotted in Fig. [2j The 
characteristic 'ribbon' shape is evident. Also shown are 
the asymptotic forms [T7] 



A(a)/£ 



f \og{£/a) + .0794, a < £■ 
n{£/a) 3 /2, a^>£. 



(5) 



There is an integrable logarithmic singularity for small a. 
For large a, the trajectories are almost circles. By com- 
bining these asymptotic limits with numerical integra- 
tion, we find the integral from Eq. ^ is J °° A 2 (a)da ~ 

2.37£ 3 . But note that J* A 2 (a)da ~ 2.3W 3 : the integral 
is completely dominated by "head-on" collisions (97% of 
the integral). 

For a sphere in three-dimensional potential flow, the 
displacement A (a) appears superficially much as for the 



cylinder in Fig. [2] but drops off more rapidly for large a: 
f | \og(e/a) - .582, a < t 



A{a)/£ 



[97r(£/a) 5 /64, 



a > I. 



(6) 



As for the cylinder, there is an integrable logarithmic 
singularity for small a, and the integral in Eq. ^ is 
aA 2 (a)da ~ .254£ 4 . The small a logarithmic singu- 
larity is mollified by the extra factor of a in the integral, 
but still J aA 2 (a)da ~ .250£ 4 : as it was for for the cylin- 
der, the integral is completely dominated by "head-on" 
collisions (98% of the integral) , due to the rapid decay of 
the displacement with impact parameter. For both the 
cylinders and spheres, the logarithmic singularity is the 
dominant contribution to the integral. The coefficient 
of the logarithm in Eqs. ^ and ^ is given by the lin- 
earized flow near the stagnation points at the front and 
rear of the cylinder or sphere, suggesting that the integral 
is easy to approximate for more complicated swimmers. 
Putting all the numerical factors together, we find the 
effective diffusivity 



1.19 Unl 3 
.266 Uni A 



(cylinders); 
(spheres). 



(7) 



We can justify these formulas dimensionally by observing 
that the frequency of collisions is linear in both U and n, 
and since we are assuming the path length is infinite the 
only other length scale is the swimmer size I. 

Direct simulation of dilute suspensions. To validate 
our theoretical predictions, we consider the encounters of 
2D swimmers moving in straight lines within the square 
S: —\L < x, y < ^L. We assume each swimmer is 
a cylinder of radius 1 = 1, with L ^> I. We initially 
place N swimmers at random positions within 5, which 
subsequently move with unit speed in a random direc- 



tion. (Figure 3(a) shows a typical initial configuration.) 



Positions are subsequently computed mod L in both di- 
rections, maintaining the number density n in S. Dilutc- 
ness requires that nl 2 < 1. A target particle initially at 
the origin moves under the potential flow created by all 
of the cylinders. Since the cylinders are typically well- 
separated, we compute their net velocity field by linear 
superposition. We show an example of this computa- 
tion in figure 3(b) for the values L = 1000, i 



io- 



The larger "ribbons" caused by drift are easily identified, 
suggesting that in this dilute limit the approximation of 
encounters as being independent will hold. 

Figure [i] shows the mean-squared displacement (|a;| 2 ) 
of a target particle over 2 x 10 6 trials (realizations) 
for N — 10 cylinders, again with unit radius and speed. 
The solid line confirms that (\x\ 2 ) grows linearly with 
time, and the dashed line shows the 2D theoretical pre- 
diction Q for cylinders. The discrepancy is due to 
Eq. |7]) only being valid in the limit of infinite dilution, 
i.e. as n£ 2 —> 0. We have verified that as the dilution 
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FIG. 3: (a) Initial configuration of N = 100 cylinders moving 
at constant speed U = 1 in random directions. Each cylinder 
has radius 1=1 and the periodic box size used in the simu- 
lation is L = 1000 (cylinders not to scale), (b) Trajectory of 
target particle, initially at the origin, integrated for 10 time 
units. 
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FIG. 4: The mean-squared displacement (solid line) of a tar- 
get particle for 2 x 10 6 realizations of ./V = 10 cylinders, with 
otherwise the same parameters as in the caption to Fig. 3(a) 



The dashed line shows the squared-displacement predicted by 
Eq. (|3f and 0, using a number density n = 10 5 . The dis- 
crepancy between the two lines is due to Q being valid in 
the limit of infinite dilution. 

is increased the slope approaches the theoretical predic- 
tion. The mean is dominated by a few trajectories with 
large displacements, corresponding to small impact pa- 
rameter a. Analogous remarks apply to a suspension of 
spheres in three dimensions: the theoretical predictions 
are verified there as well. 

Typical numerical values. We will use values for typ- 
ical krill as in [4]. We consider spheres of radius I — 
1 cm, swimming speed U = 5 cm/ sec, and number den- 
sity n = 5 x 10~ 3 cm~ 3 . Equation (JtJ) then gives an 
effective diffusivity of 7 x 10 -3 cm 2 / sec, about five times 
the thermal molecular value 1.5 x 10~ 3 cm 2 / sec, and five 
hundred times the molecular value 1.6 x 10~ 5 cm 2 / sec 
for salt. This implies a considerable enhancement to the 
molecular diffusion, but we emphasize that these values 
apply within a school of krill: the distribution and size of 
the schools themselves is a more complicated matter [7]. 
Note also that a small change in the swimmer size has a 
huge impact: for a radius of .5 cm, the effective diffusiv- 
ity is 4 x 10 -4 cm 2 / sec, an order of magnitude smaller 
than for 1 cm. If we use mean densities as discussed in 
the introduction, the effective diffusivity decreases by a 



4 



factor of 10 to 100. 

Effect of viscosity. We expect viscosity to greatly en- 
hance re. This will be the focus of future investigation, 
but for now we present a rough estimate of the impact 
of viscous no-slip boundary conditions at the surface of 
the swimmer. For inviscid flow, the displacement func- 
tion A (a) has a logarithmic singularity near the axis of 
swimming; for viscous flow near a no-slip boundary, the 
displacement function has the stronger singularity 



A (a) ~ C£ 2 /a, a -C £, 



(8) 



where C=-\/2/37rfora sphere of radius I [201 El] • This 
implies that the 3D squared-displacement integral in ([3]) 
diverges as a — > 0. The divergence of the second mo- 
ment (|£c| 2 ) is often associated with Levy flights, but here 
we are interested in scales that are much larger than the 
typical correlation length of swimming, i.e., the typical 
length A for which a swimmer travels roughly in a straight 
line before changing direction. We thus expect the over- 
all long-time transport to remain diffusive, and we can 
cap-off the displacement function at a maximum value A. 
In other words, a particle which is directly in the path of 
a swimmer cannot travel further than the swimmer itself: 
this regularizes the integral ^ to give 



2nUnt 



C£ 2 /\ 



a A 2 (a)da 



(9) 



where the lower bound of the integral is the value at 
which the displacement ^ achieves its maximum allow- 
able value, A. We introduce a transition length scale 
where we switch from the boundary-layer form A ~ a" 1 
to the inviscid form derived earlier, and find again that 
the dominant contribution to the integral arises from 
small a, as was the case for potential flow, yielding 



re~ f Un£ 4 C 2 \og(X/e). 



(10) 



For spheres, this is k ~ (2vr 3 /9)Un£ 4 log(X/£) = 
6.89Un£ 4 \og(X/£). Inserting the same numerical val- 
ues for krill as before, with a path length A = 100 cm, 
we find re ~ .8 cm 2 / sec, about 500 times the molecular 
value. Thus, including the effect of viscosity and finite 
path length has increased the effective diffusivity by a fac- 
tor of 100 over the inviscid flow case. We emphasize that 
this is a rough estimate. The path length (or swimming 
correlation length) is a measure of how much a swimmer 
tends to move in one direction before turning. Our cho- 
sen value of 1 m is not based on any evidence, but re has 
only a weak logarithmic dependence on A. Assuming a 
value A = 10 m only raises re from .8 to 1.2 cm 2 / sec. 

Any conclusion regarding the importance of biomixing 
in the oceans must be carefully qualified: at the den- 
sities inside of schools, the inclusion of viscous effects 
suggests a rather large enhanced diffusivity, comparable 
with other processes pQ, while outside of schools the ef- 
fect is much weaker. However, our viscous estimate is 



rough and more effects will need to be included to form 
a complete theory: boundary layers, more realistic shape 
distributions for the swimming bodies, wakes, spatial cor- 
relations between the swimmers, patchiness and school- 
ing, finite correlation length of swimming, distribution of 
velocities, and buoyancy and stratification effects. This 
last item is probably the most important: stratification 
can cause fluid parcels to return to their initial height af- 
ter being displaced if they can't equilibrate their density 
with their surroundings. A mechanism such as enhanced 
diffusion due to small-scale turbulence might assist this 
equilibration. 

The simplicity of our model means that prefactors and 
scalings can be computed accurately. The numerical con- 
stants we obtained depend mostly on the flow near the 
stagnation points around the swimming body. Our sim- 
ple model can serve as a platform on which to build com- 
plexity, or could be applied to other fluid-dynamical sys- 
tems where a collection of objects causes mixing, such as 
in sedimentation. 
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